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Abstract 

The mass matrix of tiiree Majorana neutrinos of equal mass is not necessarily 
proportional to the identity matrix, but expressible in terms of two angles and one 
phase. We discuss how such a mass matrix may be stable or unstable against radiative 
corrections. We then propose a model with additional explicit breaking of the threefold 
degeneracy to account for the atmospheric neutrino data, while the radiative breaking 
explains the solar neutrino data, using the large-angle Mikheyev-Smirnov-Wolfenstein 
solution. Our model requires a nonzero effective mass for neutrinoless double beta 
decay close to the present experimental upper limit of 0.2 eV. 



There are now several experimental results |T|, |^ favoring neutrino oscillations as their 
explanation. Since only the differences of the squares of neutrino masses are involved, 
there has always been a lot of theoretical interest in considering nearly degenerate neu- 
trino masses]^. The origin of their splitting may in fact be radiative]^ and some simple 
specific models have been proposed[§. In the case of 3 Majorana neutrinos, it has been 
shown [0 that all may be identical in mass and yet are inequivalent to one another because 
their mass matrix may contain 2 angles and 1 phase which cannot be redefined away. In the 
following we will discuss the stability of this construction under radiative corrections. 

We will then propose a model with additional explicit breaking of the threefold degeneracy 
to account for the atmospheric neutrino data[|l|, while the radiative breaking explains the so- 
lar neutrino datap], using the large-angle Mikheyev-Smirnov-Wolfenstein (MSW) solution[|]. 
Our model assumes all 3 neutrino masses to be of order 0.5 eV, to be consistent with their 
possible role in cosmic structure formation, in the light of recent astrophysical evidence |p.0|| 
favoring a nonzero cosmological constant. An effective v^, mass for neutrinoless double beta 



decay close to the present experimental upper limit |TlI] of 0.2 eV is also required. 



We start out with the most general 3x3 unitary matrix linking (i/g, i^fi-, ^r) with their 



mass eigenstates (i^i, z/2, z/3), i.e.[|2 

C1C3 



U 



sicse 



-iSi 
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C1C2 - SiS2S3e 



i(<53+'52— <5i) S2C3e*'^^ 
(<52-<5i) 



(1) 



-CiS2e - 5102336"^"^ C2C3 

which has 3 angles and 3 phases. Let us set S3 = (03 = 1), 5i = tt/2, 62 = S, and 
63 = 11/2 — 6, and multiply on the left by the diagonal matrix (1, —i, e"**^), then 

Cl —isi 



U 



-S1C2 
S1S2 



-IC1C2 S2e~ 
iciS2 C2e~ 



-iS 



(2) 



Since all neutrino masses are assumed equal, the mass matrix in the (z/i, 1/2, z/3) basis is 
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(3) 



just m times the identity matrix. However, in the (z/g, z/^, Ur) basis, it is given by 

Co — ■S0C2 S0S2 

M = mU*U^ = m -S0C2 -co4 + 4^^'^ S2C2(co + e^*^) 

S0S2 S2C2{co + e^'^) -Cosl + 4e^'^ 

where cq = c\ — s\ and Sq = 2s\C\. As shown in Ref. 0, this mass matrix (which has 2 angles 
and 1 phase) cannot be reduced further and there is CP violation if e^**^ 7^ =tl- Another way 
to see this is to realize that U of Eq. (2) cannot be written in the form V = DO, where D 
is a diagonal matrix containing phases and O is an orthogonal matrix. 

Even though U of Eq. (2) is nontrivial, the absence of any mass differences among the 
3 neutrinos will not result in oscillations. However, in the presence of radiative corrections 
from the charged leptons, neutrino mass differences will occur. It is easiest to do this in the 
{ui, z/2, z/3) basis. Without radiative corrections, this mass matrix is just m times the identity 
matrix. With radiative corrections coming from the wavefunction renormalization of due 



to m-r, we find 



M' 



mil 



m 



where 
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1 + 2clsle 

2csSiS2C2e -2ssCiS2C2e 
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(4) 



(5) 



with A equal to the scale at which the degenerate mass matrix is defined. 
The characteristic polynomial equation of M'/m is given by 



3 2 I r 2 2/1 22 2 2\,224 

X X ~r X\S2C2\i' ^S^l ^S^l) ' '^lCi^2 



0, 



(6) 



where x = (A — l)/2e. Hence one eigenvalue A is 1 and it corresponds to the eigenstate 



1/ 



\CsC1C2V1 - SsSiC2l'2 - S1C1S2VZ 



] = N-^/^[c2icscl - ISssDue - SiCie'^U^], (7) 



where N — Cgclcl + s^sfcl + s^c^s|. This shows that the original pattern of 2 mixing angles 
and 1 phase is generally unstable and no matter what values they take, a mass eigenstate 
exists without 1/^- as given above. For J\4' to be stable, the off-diagonal elements must go to 
zero, but since they depend on the original 2 angles and 1 phase, the condition of stability 
imposes severe constraints on their values. 

To understand the atmospheric neutrino data in terms of oscillations, we need §2 ~ ^2 ~ 
0.5, hence the requirement of stability forces si (as well as ss) to be small, which implies 
z/' ~ z/i ~ z/g. This means that we must take the small-angle MSW solution for the solar 
neutrino data. However from Eq. (6), it can be shown that the mass eigenvalues of the other 
2 states are both smaller than 1 because e of Eq. (5) is negative. Thus there cannot be any 
resonance enhancement of oscillations from interactions in the sun and this scheme does 
not work. 

We now propose to break the threefold degeneracy also explicitly by a mass term 

m'(sV2 + cVa)^, (8) 

where m' is chosen such that {m+m'y—m^ is of order 10~^ eV^ to be suitable for atmospheric 
neutrino oscillations. Then the 2x2 mass submatrix spanning ui and c'1'2 — s'v^ is given by 



M" = m 



1 -|- 2sfs|e —2s' CiSiS2C2e 

-2s'csSiS2C2e 1 + 2(c'^cfs| + s''^c^ ^s^s' dcxS2C2)e 



(9) 



Consider the simplified case of = —1 (c^ = 0), then M." is diagonal. Hence v\ and 
c!v2 — s'vj, are eigenstates with eigenvalues 

mi = m(l + 2s\s\t), ma = m[l H- 2{d'^s\c{ \ s''^cl - 2s' d S2C2Ci)e] (10) 

respectively. We choose the convention that c\ > 1/2, so that ui is mostly i/g. Therefore, 
the requirement mi < m2 means 

22^/2 22i/2 2o'/ // ' \2 

S^S2 > C S2CI + S C2 — 2s C S2C2C1 = (C S2C1 — S C2) , (11) 



resulting in the following condition on Ci: 



1 s'c'c2 + sj2si- s 

71 ^ ^ 52(1 + C'2 
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(12) 



Note that s' = is not a solution. 



With the new set of mass eigenstates, the transformation matrix U of Eq. (2) becomes 



Ci —id Si —is'si 

-S1C2 -i{c'ciC2 + s's2) -i{s'ciC2 - c's2) 
S1S2 i(c'ciS2 - s'c2) i(s'ciS2 + C'C2) 



(13) 



Using the well-known expressions for the probabilities of neutrino oscillations, we find in the 
atmospheric case, 



P^^ = l-2\U^3\\l-\U,,3\^){l-cos{Am^L/2E)) 

Pee = l-2|[/e3|'(l- |t/e3|')(l-COs(Am2L/2E)) 



fie 



= 2\Ue3\'\U^3\\l - cos{Am^L/2E)), 



(14) 
(15) 
(16) 



where Am^ = (m + m')^ — tti^- In the solar case. 



l-2\U. 



e3| 



It/esH - 2\Uei\^\Ue2\\l ' cos{A' L / 2E)) , 



;i7) 



where A'm^ = m\ — Consider again the small- angle MSW solution of the solar neutrino 
data. This requires ci ~ 1, but then Eq. (11) implies that (c's2 — s'c2)^ — |f/^3p has to be 
very small, which is disallowed by the atmospheric neutrino data. Note especially that this 



restriction is independent of the value of m. Hence only the large- angle MSW solution 113 
will be considered from now on. 

We allow e of Eq. (5) to be divided by cos^ /?, where tan (3 = V2/V1 in a two-Higgs doublet 
model. We then fix 



A/ 2 z z 

m = 1712 ~ ^1 



Am^Gprnl A 



In- 



167r2v/2cos2/5 m 



w 



22 /I / \ 2 

S1S2 - (C S2C1 - S C2) 



to be equal to 10^^ eV^ to fit the solar data. This is necessary because we will set A = 10^^ 
GeV and choose m to be 0.5 eV or less0 which then require cos^ (3 < 1 to obtain the desired 
value of A'm^. From the CHOOZ reactor data[|n], we require 



|t/e3p < 0.025. (19) 
From the atmospheric neutrino data|l[], we require 

0.84 < 4|f/^3r(l - If/MsH = sin' '^^atm < 1. (20) 
From the large-angle MSW solution||l3l to the solar neutrino datap|, we require 



0.7 < A\Uei\'^\Ue2\^ = sin' 29soi < 0.9. (21) 
From neutrinoless double beta decay |]1T|, we require 



ruee = mcQ = m{cl - sj) < 0.2 eV. (22) 

We then choose m = 0.5 eV, and scan the parameter space of s', si, and S2 for solutions. 
We find the following allowed ranges of values (each obtained with the others free): 

0.16 < s' < 0.29, 0.55 < si < 0.64, 0.66 < S2 < 0.84. (23) 

In Table 1 we show some typical solutions. It is clear that the effective neutrino mass rriee 
for neutrinoless double beta decay is not far below its current experimental upper limit of 
0.2 eV and is accessible to the next generation of such experiment s [|l5| . We note also that 



Eq. (18) depends only on the ratio m/ cos (3. Hence if we decrease m, then the same set of 
values for s', si, and S2 is a solution if we also decrease cos (3 by the same factor. Of course, 
rriee is also reduced, but if m plays a role in cosmic structure formation, then rriee cannot be 
an order of magnitude smaller than 0.1 eV. 

In the above, we have chosen = — 1 in Eq. (9). This allows Eq. (11) to be satisfied 
with the widest range of parameter values. For a general ss, Eq. (11) is replaced by 

sjsl > d'^slcl + s"c2 + 2sss' d S2C2C1 = {cs2Ci - s'c2Y + 2(1 + ss)s'c S2C2C1, (24) 
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which is clearly more restrictive. 

In conclusion, we have shown in this paper how three inequivalent mass-degenerate Ma- 
jorana neutrinos are unstable against radiative corrections due to m^. However, if we add 
an explicit mass term which breaks this threefold degeneracy to account for the atmospheric 
neutrino data, the remaining radiative splitting is able to account for the solar data, but only 
with the large-angle MSW solution, resulting in an effective neutrino mass for neutrinoless 
double beta decay close to the present experimental upper limit of 0.2 eV. 

The research of E.M. was supported in part by the U. S. Department of Energy under 
Grant No. DE-FG03-94ER40837. 
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1: Typical allowed parameter values of this model for m — 0.5 
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